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Abstract 

The thermal evolution of the hadronic parameters of charmonium in the vector channel, i.e. the J/iff resonance 
mass, coupling (leptonic decay constant), total width, and continuum threshold is analyzed in the framework of 
thermal Hilbert moment QCD sum rules. The continuum threshold Sq, as in other hadronic channels, decreases 
with increasing temperature until the PQCD threshold sq — 4 to 2 is reached at T =s 1 .22 T c {ytiq is the charm quark 
mass) and the J/iff mass is essentially constant in a wide range of temperatures. The other hadronic parameters 
behave in a very different way from those of light-light and heavy-light quark systems. The total width grows with 
temperature up to T - 1 .04 T c beyond which it decreases sharply with increasing T. The resonance coupling is 
also initially constant beginning to increase monotonically around T - T c . This behavior strongly suggests that the 
J /iff resonance might survive beyond the critical temperature for deconfinement, in agreement with lattice QCD results. 
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1. Introduction 

We discuss here the thermal evolution of the hadronic 
parameters of J /iff in the vector channel, using ther- 
mal QCD Sum Rules yj]. We refer the reader to the 
original article |01 for details. This technique has been 
used previously in the light-light and in the heavy-light 
quark sector yt]- JH, with the following emerging pic- 
ture: (i) For increasing temperature, hadronically sta- 
ble particles develop a non-zero width, and resonances 
become broader, diverging at a critical temperature in- 
terpreted as the deconfinement temperature (T c ). The 
thermal resonance broadening was first proposed in 0]. 
ii) Above the resonance region, the continuum threshold 
in hadronic spectral functions, i.e. the onset of pertur- 
bative QCD (PQCD), decreases monotonically with in- 
creasing temperature. When T — > T c hadrons disappear 
from the spectrum, (iii) This scenario is also supported 
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by the behavior of hadronic couplings, or leptonic decay 
constants, which approach zero as T — > T c . Masses, on 
the other hand, do not to provide information about de- 
confinement. 

The thermal behavior of the heavy-heavy quark corre- 
lator should be different from that involving at least 
one light quark since: a) In the light-light and heavy- 
light quark sector, the PQCD contribution is domi- 
nated by the time-like spectral function (annihilation 
term), which is relatively unimportant in relation to the 
light quark condensate contribution, being the scatter- 
ing PQCD spectral function highly suppressed. Instead, 
for heavy-heavy quark systems this term becomes in- 
creasingly important with increasing temperature while 
the annihilation term only contributes near threshold; 
b) The non-perturbative QCD sector in the operator 
product expansion (OPE) of light-light and heavy-light 
quark correlators is driven by the light quark conden- 
sate, responsible for the behavior of the continuum 
threshold since s (T)/s (0) * {{qq))/{qq) fll-fl. The 
light quark condensate is the order parameter for chi- 
ral symmetry restoration. In contrast, for heavy-heavy 
quark correlators the leading power correction in the 
OPE is that of the gluon condensate, which has a very 
different thermal behavior. In this approach, the criti- 
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cal temperature for deconfinement is a phenomenologi- 
cal parameter which does not need to coincide with e.g. 
the critical temperature obtained in lattice QCD |8fl. In 
fact, results from QCD sum rules lead to values of T c 
somewhat lower than those from lattice QCD. In order 
to compare with other approaches, we express our re- 
sults in terms of the ratio T/T c . 

We find for charmonium in the vector channel that the 
continuum threshold, Sq(T), decreases with increasing 
T, being driven by the gluon condensate and the PQCD 
spectral function in the space-like region, until it reaches 
the PQCD threshold s = 4m 2 Q at T - 1.22 T c (m Q is 
the charm quark mass). Below this value of sq the sum 
rules cease to be valid. The J/if/ mass remains basically 
constant as in the light-light or heavy-light systems. We 
have, however, a very different thermal evolution of the 
width and the coupling. Both are almost independent of 
T up to T - 0.8 T c where the width begins to increase 
substantially, but then above T - 1 .04 T c it starts to de- 
crease sharply, and the coupling increases also sharply. 
This suggests the survival of the J/iff resonance above 
the deconfinement temperature. 

The PQCD spectral function in the space-like region 
plays here a very important role. Non-relativistic ap- 
proaches to charmonium at finite temperature would 
normally miss this contribution. In fact, the complex en- 
ergy plane in the non-relativistic case would only have 
one cut along the positive real axis, which would cor- 
respond to the time-like (annihilation) region of PQCD. 
The space-like contribution (q 2 = {or - |q| 2 ) < 0) in the 
form of a cut in the energy plane centered at the origin 
for -|q] < lo < |q|, would not be present in the non- 
relativistic case. 



2. Hilbert Moment QCD Sum Rules 

We consider the correlator of the heavy-heavy quark 
vector current at finite temperature 

IV(<7 2 , T) = -(g flv q 2 - q M q v ) Il{q 2 , T) 

= ij d*Xe** 9(x ) « |[V„(jc) , V?(0)]| », (1) 

where V M (x) =: Q(x)y t ,Q(x) :, and Q(x) is the heavy 
(charm) quark field. The matrix element above is under- 
stood to be the Gibbs average in the quark-gluon basis. 
The imaginary part of the vector correlator in PQCD 
at finite temperature involves two pieces, one in the 
time-like region (q 2 > 4wg), Im Tl a (q 2 , T), which sur- 
vives at T=0, and one in the space-like region (q 2 < 0), 
Im Tl s (q 2 , T), which vanishes at T=0. To leading order 
in PQCD we find 



-ImIl a (q 2 ,T) = -J-j 
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where v = 1 - Am Q /q , wiq is the heavy quark mass, 
q 2 — cj 2 - q 2 > 4wg, and np(z) = (1 + e z )~ x is the Fermi 
thermal function. In the rest frame of the thermal bath, 
|q| — > 0, the above result reduces to 

1 v (3 _ v 2) 

- Imll a (a>, T) = \ - [1 - 2n F (to/2T)] 

n in 1 

X 0(u - 2m Q ) . (3) 

The quark mass is assumed independent of T, a good 
approximation for T < 200 MeV |@1 • Only the leading 
order in the strong coupling will be considered here. 

The PQCD piece in the space-like region demands a 
careful analysis. In the complex energy plane, and in the 
space-like region, the correlator FI(g 2 ), Eq.(l), has a cut 
centered at the origin and extending between to = -|q| 
and a> — |q|. In the rest frame this cut produces a delta 
function 5(a> 2 ) in the imaginary part of fl(<7 2 ). The result 
is 

1 2 , , 
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The corresponding hadronic representation is 
parametrized in terms of the ground state resonance,the 
J '/iff, followed by a continuum given by PQCD after a 
threshold so > M 2 . In the zero width approximation, 
the hadronic spectral function is 

- ItnUis, T)\ HAD = -ImTl(s, T)\ RES 6(s - s) 
n n 

+-ImU(s, T)\ PQCD 6(s - s ) 

7T 

= 2f 2 (T)6(s - M 2 V (T)) + -ImUis, T) a 0(s - s ) , 
n 

(5) 

where s = q 2 = at 2 - q 2 . The leptonic decay constant is 
defined as < 0|V M (0)|V(£) >= V2 M v f v e M . 

When considering a finite (total) width the following 
replacement will be understood 
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The hadronic scattering term, due to current scatter- 
ing off D-mesons, is negligible J2l- The correlation 
function Tl(q 2 ,T), Eq.(l), satisfies a once subtracted 
dispersion relation. To eliminate the subtraction one can 
use Hilbert moments, i.e. 



9 (-) / d \N , 



i r°° 

~ n Jo (s + 



ds 



2VV+1 



Q 2 ) 



Im Ii{s,T) , 



(7) 



where N = 1,2,..., and Q 2 > is an external four- 
momentum squared, to be considered as a free parame- 
ter. Using Cauchy's theorem in the complex s-plane, the 
Hilbert moments become Finite Energy QCD sum rules 
(FESR), i.e. 



MQ 2 , T)\ RES = ip N (Q\ T)\ QCD , 



(8) 



where 
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Im Tl a (s, T) 



Im n s (s, T) + ip N (Q 2 , T)\ NP , 
(10) 

and tallfs, T)\res is given by the first term in Eq.(5) 
modified in finite-width according to Eq.(6), and the 
PQCD spectral functions are given by Eqs.(3) and (4). 

The dimension d=4 non perturbative term in the OPE 
is well known in the literature, see yfl for details. The 
dependence on N is quite cumbersome and it is propor- 
tional to the gluon condensate W^G 2 \\. At low tem- 
peratures, this condensate 

has been calculated in chiral perturbation theory fioll . 
In this framework the condensate remains essentially 
constant up T ~ T c - 100 MeV, after which it de- 
creases sharply. In order to go beyond the low temper- 
ature regime of chiral perturbation theory, lattice QCD 
provides the right tool. A good approximation ll ill is 
given by the expression 
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where T* w 150 MeV is the breakpoint tempera- 
ture where the condensate begins to decrease apprecia- 
bly, and T* c * 250 MeV is the temperature at which 

«fG 2 »r c = °- 

Returning to the Q 2 dependence of the Hilbert moments, 
Eq.(7), we shall fix Q 2 and sq{0) from the experimental 
values of the mass, the coupling, and the width at T=0. 
At finite temperature there are non-diagonal (Lorentz 
non-invariant) condensates that might contribute to the 
OPE. Non-gluonic operators are highly suppressed 01, 
ni2fl so that they can be safely ignored. We have con- 
sidered also a gluonic twist-two term in the OPE intro- 
duced in lfl3ll . and computed on the lattice in fl4ll . Its 
impact is small,(2-6)%, and plays no appreciable role in 
the results. 

3. Results 

We begin by determining Sq and Q 2 at T=0 from the 
moments, Eq.(8), and using as input the experimental 
values |H1 My = 3.097 GeV, f v = 196 MeV, and 
IV = 93.2 keV, as well as m Q = 1.3 GeV, and H3] 
(0|^G 2 |0> - 5 x 10~ 3 GeV 4 . In the zero-width approx- 
imation one finds from Eq.(9) that 



>fl(Q 2 )\RES _ fl(Q 2 )\RES 
n(Q 2 )\RES >fi(Q 2 )\RES 



(12) 



Given the extremely small total width of the J/ijj it 
turns out that the above relation also holds with extreme 
accuracy in finite width. Using Eq.(8) this leads to 



ifl(Q 2 )\QCD _ <P2(Q 2 )\qCD 
*P2(Q 2 )\qCD >P3(Q 2 )\qCD ' 



(13) 



which depends only on the two unknowns so and Q 2 , 
and provides the first equation to determine this pair of 
parameters. The second equation can be e.g. Eq.(8) 
with N = 1. In this way we find that Sq = 11 .64 GeV 2 , 
and Q 2 = 10 GeV 2 reproduce the experimental values 
of the mass, coupling, and width of J/tfr within less than 
1%. This whole set of hadronic parameters will then be 
used to normalize the corresponding parameters at finite 
temperature. In this way, see 12J] for details, we were 
able to find the thermal evolution of So, the J/ifi mass, its 
width and its coupling. We show here only the behavior 
of the width and the coupling (Figs. 1 and 2) since these 
are the most important results of this analysis. 

Both the width and the coupling can only be deter- 
mined up to Tf =! 1.1 T c beyond which so(T) < My(T) 
and the FESR integrals have no longer a support. The 
temperature behavior of the width and the coupling 
shown in Figs. 1 and 2 strongly suggests the survival 
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Figure 1: The ratio rV(7")/rV(0) as a function of T/T c . 



Figure 2: The ratio fv(T)l /v(0) as a function of T/T c . 



of the Jj\p above the critical temperature for decon- 
finement. This conclusion agrees with results from 
lattice QCD H, but disagrees with non-relativistic 
determinations. As pointed out earlier, the reason for 
this disagreement might very well be the absence of the 
central cut (QCD scattering term) in the energy plane 
in non-relativistic frameworks. 
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